Role of low-/ component in deformed wave functions near the continuum threshold 
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The structure of deformed single-particle wave functions in the vicinity of zero energy limit is 
studied using a schematic model with a quadrupole deformed finite square-well potential. For this 
purpose, we expand the single-particle wave functions in multipoles and seek for the bound state 
and the Gamow resonance solutions. We find that, for the K n — + states, where K is the z- 
component of the orbital angular momentum, the probability of each multipole components in the 
deformed wave function is connected between the negative energy and the positive energy regions 
asymptotically, although it has a discontinuity around the threshold. This implies that the K n = + 
resonant level exists physically unless the I = component is inherently large when extrapolated 
to the well bound region. The dependence of the multipole components on deformation is also 
discussed. 

PACS numbers: 21.10.Pc, 21.60.Jz, 24.30.Gd 



I. INTRODUCTION 

Physics of nuclei located far from the (3 stability line 
has been one of the main current subjects of nuclear 
physics. One of the unique properties of drip- line nuclei 
is that the Fermi level lies close to zero. Understand- 
ing of single-particle levels in the continuum is essential 
in describing the nuclear structure close to, and beyond, 
the drip line, since the shell structure of both bound and 
continuum levels plays an important role in many-body 
correlations such as deformation and pairing. 

It has been argued recently that, as the binding en- 
ergy approaches zero, the s-wave component of a bound 
single-particle wave function behaves uniquely in a de- 
formed potential, and plays a dominant role in Nilsson 
levels with Q n = l/2 + [l],|3- Naively, resonant levels can 
be considered as an extension of bound states into the 
positive energy regime. Therefore, if the s-wave com- 
ponent keeps dominant in the continuum, the level with 
SV = 1/2+ might not exist as aphysical state. Notice 
that, for a Nilsson HamiltonianQ, single-particle levels 
with f2 = 1/2 belonging to high-j orbit comes down in en- 
ergy in a prolately deformed potential. These states play 
an important role in generating the deformed shell struc- 
ture. It is therefore crucially important to investigate the 
role of \ow-l component in a deformed wave function for 
£1 = 1/2 states and its transition from bound to resonant 
levels. 

The structure of deformed single-particle levels in the 
continuum has been investigated in a few publications. In 
Ref. 4] , the resonance energy of negative parity states was 
studied by employing the Gamow wave function. The 
Analytic Continuation in the Coupling Constant (ACCC) 
method was applied to study single-particle resonance 
states in spherical and deformed nuclei Jjj. Using the 
multi-channel scattering approach, Ref. |(j has studied 
how the single-particle energies change from bound to 
resonant levels when the depth of the potential is var- 
ied. In order to fully understand the structure of de- 



formed single-particle levels in the continuum, however, 
a detailed study of the wave function components is still 
necessary, in addition to the resonance energy itself. 

In this paper, we investigate the structure of deformed 
wave functions around zero energy using the Gamow 
state representation for a resonant state. To this end, 
we use a schematic model: a Y20 deformed finite square- 
well potential without spin-orbit force. This enables us to 
determine the single-particle wave function analytically. 
To use the Gamow state for resonance has a certain ad- 
vantage in analyzing the deformed wave function. That 
is, we are able to treat the bound and the resonant lev- 
els on the same footing, because the Gamow states are 
normalizable just like the bound states 0. It is then 
straightforward to see how the fraction of each compo- 
nent in the deformed wave functions changes when the 
single-particle level changes its character from bound to 
resonant. A slight disadvantage of this approach is that 
the expectation value with the Gamow states, including 
the probability of wave function components, becomes 
complex numbers. However, this is not a big defect for 
our purpose, since the physical quantity of the expecta- 
tion values can be obtained by taking their real part|^,@. 

The paper is organized as follows. In the next sec- 
tion, we present our model for a deformed single-particle 
potential. Numerical results and discussion are given in 
Sec IIIII Finally, we summarize the paper in Sec II VI 



II. MODEL 

Our purpose is to study the structure of wave function 
in a deformed single-particle potential. To this end, we 
employ a schematic model for the single-particle poten- 
tial, that is, a deformed square- well potential without the 
spin-orbit force, 
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where R(r) — Ro(l + ^2^o(^))- For simplicity, we ex- 
pand this potential up to the first order of deformation 
parameter fa and obtain 

V(r) ~ -Vo [0(R - r) + R a l3 2 Y 20 (f)5{r - R Q )] . (2) 

In order to solve the Schrodinger equation with this po- 
tential, we expand the wave function in the multipoles as 



tt*:(r) = J2 



ui K (r) 



(3) 



where the quantum number K(= A) is the ^-component 
of the orbital angular momentum I. By projecting out 
each multipole component, we obtain the coupled equa- 
tions for the radial wave functions given by 
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For the positive energy solution, E > 0, we impose the 
boundary condition corresponding to the Gamow state 
for resonance. That is, the wave function is regular at 
the origin and satisfies the out-going boundary condition 
u(r) ~ e lkr asymptotically. This boundary condition is 
satisfied only if the energy is complex, E — h 2 k 2 /2m = 
Eji — iT/2, where En and T are the resonance energy and 
the width, respectively. In the case for T = and Er < 
0, the Gamow state wave function is equivalent to the 
bound state wave function, which satisfies the decaying 
asymptotics u(r) ~ e~ 7r , where 7 = \J ' —2mE^/h 2 . 

The solutions of the coupled-channels equations @ 
therefore read (we omit the subscript K for simplicity 
of notation) , 
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where k x = v / 2m(E + V )/h 2 , k = ^/2mE/h 2 , and 

ji(x), h\ + \x) are the spherical Bessel and Hankel func- 
tions, respectively. The amplitudes Ai and Bi are deter- 
mined by the matching condition at r = R given by, 
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ui(R-) =ui(R+), 
[tt{(Af)-«Kfl-)] =V o Ro02 
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where R± represents lim e ^o -Ro ± £■ 

The bound state wave function is normalized as 



1= j dr\V K (r)\ 2 = J2 N i> 



where 



N, = 



dr\ui(r)\ z 



(8) 



(9) 



1.2 £,=0.5, i? =5.0(fm) 



> 

CD 



1 

0.8 
0.6 
0.4 
0.2 




K n =0 + 



-4-3-2-10123 
E R (MeV) 

FIG. 1: The real part of the energy and the resonance width 
for a K w = + state with various potential depths. The 
corresponding potential depths are shown in Fig|5| 



The Gamow state wave function can be also normal- 
ized by introducing the regularization factor as Zel'dovich 
proposed^] 



Ni — lim / dre 

e->0 
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Using a property of the spherical Bessel function [Til], one 
can evaluate the first term as 



driAirjiikxr)} 2 
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{ji(hRo)} 2 - jt-i(kiRQ)ji +l (kiRo)y (12) 



The second term can be also evaluated using the con- 
tour integral method or equivalently the Complex Scaling 
Method (CSM). The result is given by [l^. 
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z~ er2 {Bi rhl +) (kr)} 2 
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Note that the fraction of multipole components Ni is in 
general a complex number for the Gamow state wave 
function. 



III. RESULTS AND DISCUSSION 

Let us now discuss the behaviour of the low-i com- 
ponents in deformed wave functions. In Sec. IIII AI we 
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FIG. 2: The real part of the energy for JT = 0+ state as a Ener 9y ( MeV > Ener 9y < MeV ) 

function of the potential depth. In the inset, the behaviour 
around zero energy is enlarged. 

FIG. 4: The eigenphase sum and its energy derivative for 
K n = 0+ state with potential depth V = 39.0 MeV (the left 
panel), and Vb = 41.6 MeV (the right panel). 
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FIG. 3: Same as Fig0 but in the logarithmic scale. The 
solid line is an expectation for the pure s-wave configuration 
given by Eq. flU . 



vary the potential depth for a fixed deformation parame- 
ter, while we vary the deformation parameter for a fixed 
potential depth in Sec. IIII Bl 



A. Dependence on potential depth 

We first study the wave functions at a fixed deforma- 
tion, /?2 = 0.5. Figure ^ shows the real and imaginary 
parts of the energy for a K w = + state in varying the 
potential depth Vq . The correspondence between the po- 
tential depth and the real part of the energy is shown 
in Fig[3 We observe that the width is quite large even 
for a small values of positive energy. This large width is 
caused by the admixture of the I = component in the 
wave function. Indeed, as shown in Fig|3 in the small 
positive energy region (0.1 MeV < U(E) < 1.0 MeV), 
the behavior of the width is consistent with the relation 
expected for the s-wave resonance state 0, S > 



T oc U(E) l+1 ^ x 3?(7V ; ; 



2=0 



(14) 



where Sft(-E) denotes the real part of E. 

Below En —0.1 MeV, the width is larger than the solid 
line, which predicts r = at Er — 0. Also, we did not 
find a physical solution between Vq= 41.62 and 41.68 
MeV, as is shown in the inset of FigJ2J These might be 
related to the possible presence of the anti-bound and 
'crazy' resonance states, as presented in Ref. ^ij for a 
spherical square- well potential (see Fig.l of Ref. [t3|). 
Above 1.0 MeV also, the width is larger than that ex- 
pected by Eq.jnj. This is due to the fact that the rela- 
tion Ea. (|14fl is valid only for small values of k |l3j . 

In Fig|3 we see that the slope of the single-particle 
energy as a function of the potential depth, dE/dVo, or 
equivalently dE/dA, where A is the mass number, be- 
comes smaller in approaching the zero binding energy. 
For a spherical square-well potential, it has been shown 
that dEi/dA — > for I = in the limit of zero bind- 
ing [l^ . This is due to the fact that the s-wave function 
can be easily extended outside the nuclear potential and 
also the kinetic energy is reduced due to the absence of 
the centrifugal barrier ^^|. This property implies that 
the / = component becomes dominant in a deformed 
wave function around the zero-binding region. On the 
other hand, the slope has a finite value in the positive 
energy region even in the limit of zero energy, thus the 
slope has a discontinuity around zero energy. Therefore, 
a care must be taken, as discussed in Ref. [l4j , when one 
estimates the energy of a deformed resonant level with 
K* = 0+ by using the ACCC method 0. 

The resonance energy and width can be also estimated 
using the eigenphase sum A(E) [l^. It is defined in terms 
of the eigenvalues of the scattering matrix (S'-matrix) as 

(tfSU) aa , = e 2i5 «W5 aia ,, A(E) = J2$a(E). (15) 

a 

The resonance energy and width are identified with 
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FIG. 5: The real part of the fraction for each multipole com- 
ponent Ni for the K n = + state. The solid, dotted, and 
dot-dashed lines indicate the 1 — 0,2 and I — 4 components, 
respectively. 
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FIG. 6: Same as Fig[l] but for a K n = 0~ state. 



the peak energy of dA(E)/dE and its FWHM, 
respectively 17J. Figure 0] shows the eigenphase sum for 
the K n = + state with two different potential depths. 
Comparing Figsni2land0] we see a good correspondence 
between the two definitions of resonance state, i.e., the 
Gamow state representation and the approach with the 
eigenphase sum. 

We now discuss the energy dependence of the fraction 
of the multipole components in the deformed wave func- 
tion. Figure [5] shows the real part of the fraction for each 
multipole component in the Gamow state wave function 
with K v = + . When the binding energy approaches 
zero, the s-wave component in the deformed wave func- 
tion becomes dominant. In contrast, in the positive en- 
ergy region, all the multipole components have a finite 
value even in the zero energy limit and show similarity 
with the well bound cases. As we will discuss in the next 
section (see FigllOlbelowh the state shown in Fig[5]origi- 
nates from the 2d orbit in the spherical limit. This states 
couples with the lower- lying 2s, lg and the higher- lying 
3s states. The dominant component is I = 4 both in well 
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FIG. 7: Same as FigHJ but for the K* = 0" state. 



bound and in resonant levels, as one sees in FigOO This 
suggests that both the well bound and the resonant levels 
have a similar property to each other and the intuitive 
picture that the resonant level is an extension of a bound 
level into the continuum is valid. 

Only at the limit of zero binding, the singular behav- 
ior of the I = component appears. This is entirely due 
to the property of the normalization integral, Ea. 1)13(1 . 
Since the Gamow state wave function is equivalent to 
the bound state wave function for En < 0, V = 0, Ea. l|13|) 
holds both for the resonance and the bound states. For 
small values of k, Ea. (ll3fl is proportional to k 21 " 1 as dis- 
cussed in Refs. pi UJ$, that diverges only for I = as 
k — > 0. When the total wave function ^ k is normalized 
according to Eq.JSJ), then only the s-wave component is 
allowed in the wave function 0. This condition is al- 
ways met for the bound state when the binding energy 
approaches the threshold. In principle, the same consid- 
eration can apply also to the resonance state when the 
resonance energy approaches zero from the positive en- 
ergy side. However, as we show in Fig^ the resonance 
state acquires a relatively large width even when the real 
part of the energy is infinitcsimally small. Since k is de- 
fined as k = y / 2m(E R - iT/2)/h 2 , it remains a constant 
even if En itself approaches zero. This leads to the dis- 
appearance of the "s-wave dominance" in the positive 
energy side. 

We next study the case for K v = 0~. In Figs. Elandd 
we show the dependence of the single-particle energy on 
the potential depth. In contrast to the case for = + , 
due to the presence of the centrifugal barrier, we do not 
see any singular behavior around zero energy. Single- 
particle energies are connected smoothly when changing 
the potential depth, and the width increases gradually 
in the small positive energy region. Figure |S] shows the 
fraction of each multipole component in the Gamow state 
wave function. As the binding energy approaches zero, 
the p-wave component becomes relatively large, that is 
consistent with the dominance of low-i component in the 
limit of zero binding energy discussed in Ref. 2]. The 
fractions are connected smoothly and asymptotically in 
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FIG. 8: The real part of the fraction for each multipole com- 
ponents Ni for the K n — 0~ state. The solid, dotted, and 
dot-dashed lines indicate the 1 — 1,3 and I — 5 components, 
respectively. 
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FIG. 10: Single-particle energies for the K w = + state as a 
function of deformation parameter p2- The potential depth is 
Vo = 45.0 (MeV), and the potential radius R = 5.0 (fm). 
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FIG. 9: Same as FigQ except for varying the deformation. 
This level corresponds to the one originating from the 2d orbit 
in the spherical limit. The deformation dependence of the 
single-particle energies is shown in Fie llOl 



the bound and resonant regions. 



B. Deformation dependence 

In this subsection, we study the deformation depen- 
dence of the low-Z component in deformed wave functions 
for a fixed potential depth. In the realistic situation, the 
location of single-particle levels changes as a function of 
nuclear deformation. Especially, the levels of = 1/2 
(K = 0) with (without) spin-orbit force belonging to 
high-j (high-/) orbit in the spherical limit play an im- 
portant role in nuclear deformation. 

Figure shows the resonance energy and width when 
the deformation parameter is varied from @2 = 0.0 to 0.5. 
The potential depth Vq and the radius Rq are set to be 
45.0 MeV and 5.0 fm, respectively. This state belongs to 
the 2d orbit at fa = 0.0 as shown in Fig^| At around 
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FIG. 11: Same as Figl5*l except for varying the deformation. 



zero energy, we see the similar behavior as in Fig^ the 
width is quite large even for the small values of posi- 
tive energy, which implies that the I — component is 
responsible for the width of the resonant level. 

The corresponding wave function components for this 
state are shown in Fig^] As in Fig|Sl we see the 
singular behavior for the s-wave component at around 
zero-binding energy, corresponding to the "s-wave dom- 
inance" in the limit of zero binding. Except for the 
zero-energy region, however, we see that the fraction of 
each multipole components is linked asymptotically and 
is smoothly connected to the d-state resonant level in 
the spherical limit. From this calculation, it is evident 
that the singular behavior of the I = component for 
the K 71 = + state occurs only just below the continuum 
threshold and this state is connected to the physical reso- 
nant level in the continuum. Furthermore, the fraction of 
each-/ components in the deformed wave function is con- 
nected smoothly from the bound to the resonant levels 
except for the region near the threshold. 
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IV. SUMMARY 

We have analyzed the structure of the deformed wave 
functions around zero energy using the Gamow state 
wave function for resonance, with which one can treat the 
resonant and bound levels on the same footing and thus 
analyze the wave function continuously from the nega- 
tive to the positive energy regions. For this purpose, we 
developed a schematic model with a deformed square- 
well potential. Since the wave functions can be obtained 
analytically with this model, detailed analyses of the de- 
formed wave functions were possible. For a K 1 * = + 
state, we have found a singularity in the resonance width 
as well as in the s-wave component in the deformed wave 
function at around zero energy. That is, the width be- 
comes considerably large even in the small positive en- 
ergy region and the I = component approaches unity 
in the limit of zero binding. We have shown that the 
"s— wave dominance" occurs only at the threshold of con- 
tinuum. Far from the zero energy region, the probability 
of each-^ components is connected asymptotically. This 
implies that the K" = + resonant level exists unless the 
I = component is large inherently when extrapolated 



to the well bound region. In contrast, for the K n = 0~ 
state, we did not find any singular behavior even in the 
zero-energy limit. The single-particle energies are con- 
nected smoothly when changing the potential depth, and 
the width increases gradually in the small positive energy 
region. The probability of each-/ component in the wave 
function is also connected smoothly and asymptotically 
between the bound and the resonant regions. 
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